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ABSTRACT 

We critically examine the dependence of spectral index on luminosity in optically se- 
lected AGN samples. An analysis of optically selected high-redshift quasars showed 
an anti-correlation of aox, the spectral index between the rest-frame 2500 A and 2 
keV, with optical luminosity (Miyaji et al. 2006). We examine this relationship by 
means of Monte Carlo simulations and conclude that a constant aox independent of 
optical luminosity is still consistent with this high-z sample. We further find that that 
contributions of large dispersions and narrow range of optical luminosity are most im- 
portant for the apparent, yet artificial, aox —lo correlation reported. We also examine 
another, but more complete low-z optical selected AGN sub-sample from Steffen et 
al. (2006), and our analysis shows that a constant aox independent of optical lumi- 
nosity is also consistent with the data. By comparing X-ray and optical luminosity 
functions, we find that a luminosity independent aox is in fact more preferred than 
the luminosity dependent aox model. We also discuss the selection effects caused by 
flux limits, which might systematically bias the — lo relation and cause discrepancy 
in optically selected and X-ray selected AGN samples. To correctly establish a de- 
pendence of aox of AGNs on their luminosity, a larger and more complete sample is 
needed and consequences of luminosity dispersions and selection effects in flux limited 
samples must be taken into account properly. 

Key words: galaxies: active ~ quasars: general - X-rays: galaxies - methods: statis- 
tical. 



1 INTRODUCTION 

The dependence of the spectral index aox of active galac- 
tic nuclei (AGN) on redshift and luminosity has important 
astrophysical implications on AGN evolution and thus has 
been studied for many years (e.g. Avni & Tananbaum 1982 
Wilkes et al. 1994; Green et al. 1995; Bechtold et al. 2003 
Vignali et al. 2003a; Strateva et al. 2005; Steffen et al. 2006 
Hopkins et al. 2007; Kelly et al. 2007). aox is defined as 



aox 



log(/2 kov//^, 
l0g(i^2 kev/2^„ 



(1) 



where f2 keV and f = are the rest-frame flux densities at 

Ktv •'2500 A 

2 keV and 2500 A, respectively. Dependence of aox on red- 
shift means evolution of the accretion process in cosmic time. 
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Most studies have concluded that there is no evidence for 
a dependence of aox on redshift (e.g. Avni & Tananbaum 
1982; Strateva et al. 2005), although some studies found that 
aox is correlated with redshift (Bechtold et al. 2003; Kel- 
ley et al. 2007). Dependence of aox on luminosity means 
a non-linear relationship between X-ray and optical lumi- 
nosity {Lx oc Lo, e 7^ 1), which provides insight into the 
radiation mechanism. An anti-correlation between aox and 
the optical luminosity has been found by many authors in 
optically selected AGNs with follow-up X-ray measurements 
at some different epoch, which means that these AGNs span 
a larger range in optical luminosity than in X-ray luminos- 
ity (e.g. Vignali et al. 2003a; Strateva et al. 2005; Miyaji et 
al. 2006, hereafter M06; Steffen et al. 2006, hereafter S06). 
Meanwhile, whether aox depends on luminosity in X-ray 
selected AGN samples remains unknown (Hasinger 2004; 
Frank et al. 2007). 

As pointed out by Yuan et al. (1998), one of the prob- 
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lems in such studies is that an apparent, yet artificial corre- 
lation between aox and optical luminosity can be caused by 
dispersions in the optical luminosity. In section 2, we present 
analysis of the relationship between aox and optical/X-ray 
luminosity using data presented in M06, in which wo find 
that dispersions in luminosity can be entirely responsible 
for the claimed dependence. We also discuss the determin- 
ing factor for this behavior. 

Another problem is the degeneracy between redshift 
and luminosity in flux-liiiiitcd samples, where redshift and 
luminosity are strongly correlated. In section 3, we examine 
a sub-sample from S06 containing 187 AGNs, which more 
completely fills the redshift and optical luminosity plane and 
thus is less affected by such degeneracy. In section 4, we com- 
pare the optical quasar luminosity function from Richards 
et al. (2006) with X-ray quasar luminosity function from 
Barger et al. (2005) in different aox models. In section 5, 
wo discuss the selection effects in flux limited samples and 
the consequently discrepancy in optically selected and X-ray 
selected samples. Discussion and conclusions are presented 
in section 6. 

We mostly use the logarithms of luminosities and de- 
note them as Ix = log L2 kcv and lo = log I/2500 , then 
aox = 0.3838(ia, — lo), where L2 kev is the 2 keV monochro- 
matic luminosity and L2500 is the 2500 A monochromatic 
luminosity in units of erg s~^ Hz~^. We adopt the cur- 
rently favored cosmology model with Hq=7Q km s"'^ Mpc^^, 
Q.M = 0.3, and = 0.7 (e.g. Spergel et al. 2007). 
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Figure 1. 61 high-redshift quasars in Miyaji et al. (2006). Cir- 
cles indicate X-ray detected quasars, while arrows indicate up- 
per limits. Panel (a) olqx dependence on the rest frame 2500 A 
monochromatic luminosity. The solid line indicates linear regres- 
sion results from the EM Algorithm in ASURV. Panel (b) rest 
frame 2 keV monochromatic luminosity against 2500 A one. The 
solid line indicates the ILS bisector result, and dashed lines in- 
dicate 1LS(Y|X) and 1LS(X|Y) results from the EM Algorithm 
in ASURV, respectively. The e = 1 relation is shown by dotted 
line for comparison. Panel (c) aox dependence on the rest frame 
2 keV monochromatic luminosity. The solid line indicates linear 
regression result using X-ray detected quasars. Spearman corre- 
lation coefficients and slopes of fitting lines are indicated in the 
upper right panel. 



2 ANALYSIS OF MIYAJI ET AL. (2006) 
SAMPLE 

2.1 Data Analysis 

The sample we use in this section is consisted of 61 high- 
redshift {z > 2.9) quasars in Figure 2 of M06 (Miyaji et al. 
2006; Vignali et al. 2003b, 2005), excluding the three quasars 
from archival Chandra data in Table 3 in M06 which might 
be biased toward higher X-ray fluxes. Only six of thcni have 
no X-ray detection. M06 found a correlation of aox with 
optical luminosity for this high-z sample, while they kept 
the discussion on the lo — aox relation open because of 
possible optical selection effects for variable AGNs which 
preferentially pick up the optically brighter phases. 

To illustrate how dispersions produce artificial correla- 
tions in this sample, we carry out two independent analysis. 
The flrst one is linear regression of aox, lo and Ix in ob- 
served data. Without further description, we perform linear 
regression using methods as follows throughout the paper: 

(i) For the aox — lo correlation, we use the EM algorithm 
in ASURV (Isobe, Feigelson, & Nelson 1986) to derive linear 
regression parameters, including X-ray undetected quasars; 

(ii) The EM linear regression algorithm in ASURV is 
based on the traditional ordinary least-squares method 
which minimizes the residuals of the dependent variable 
(OLS(Y|X)). However, for l^ — lo correlation, both vari- 
ables are observed and a different result can be obtained if 
residuals of the independent variable arc instead minimized 
(e.g. S06). Following S06, we perform linear regression with 
ASURV using EM algorithm, treating as the dependent 
variables (ILS(Y|X)) and treating to as the dependent vari- 



ables (ILS(X|Y)), then use the equations given by Isobe et 
al. (1990) to calculate the bisector of the two regression lines. 

(iii) For the aox — Ix, where both independent and de- 
pendent variables arc upper limits, only Schmitt's binned 
method in ASURV is available which may suffer from sev- 
eral drawbacks (Sadler et al. 1989). Hence we abandon upper 
limit points in aox — Ix plane and only use X-ray detected 
quasars to derive linear regression parameters. 

Spearman correlation coefficients are calculated using 
ASURV including X-ray undetected quasars. Observational 
data together with linear regression slopes and Spearman 
correlation coefficients for aox — lo, Ix — h and aox — Ix 
correlations are shown in Figure 1. Conflicting correlations 
arise due to dispersions in luminosity: aox = — O.lS/o + 
const as shown in solid line in Panel (a), but aox = 0.20lx + 
const as shown in solid line Panel (c). Therefore the same 
data produce two totally different results: Ix = 0.53lo+const 
or Ix = 2.08^0 + const, i.e. e < 1 or e > 1 if Lx oc Lq. 
Therefore, depending upon how the regression is done, the 
conclusion on the relationship between the optical and X- 
ray luminosities can be significantly different. As shown in 
Panel (b), the slope of the lx — lo relation does depend on 
which luminosity is used as the dependent variable. When 
treating lx as the dependent variable, the slope is 0.54±0.14 
(the flatter dashed line), while it changes dramatically to 
2.08 ± 0.45 (the steeper dashed line) when treating lo as the 
dependent variable. Using ILS bisector, we find the slope to 
be 1.05 ± 0.20 (solid line), which is consistent with e = 1 
(dotted line). 

The second part of analysis is done with Monte Carlo 
simulations. Following Yuan et al. (1998), we assume intrin- 
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sic optical and X-ray luminosities lo and Ix with a constant 
mean aox, 

L = lo + aox /0.3838, (2) 

and dox = —1.728 is given by the mean of observed 
0.3838(Zx — lo) using the Kaplan- Meier estimator in ASURV, 
including the six quasars with upper limits. The above re- 
lationship is plotted as dotted line in Panel (b) of Figure 1. 
The observed optical and X-ray luminosities are assumed to 
be the intrinsic luminosities modified by independent Gaus- 
sian dispersions 

lo = lo + Slo, lx^L+ SL = L + aox/0.3838 + 51^, (3) 

where 5lo and Six are Gaussian distributed dispersions with 
standard deviations Uo and Ox respectively. Thus, the dis- 
tribution of aox is Gaussian with standard deviation 

<T„o^ =0.3838(^T^-hcT^)'''^ (4) 

where Oaox ~ 0.154 is the standard deviation around the 
linear relationship of Equation (2) for this sample, using 55 
X-ray detected AGNs. The ratio of the standard deviations 
of the optical to the X-ray luminosity dispersion is defined 
as 

R. = ^. (5) 

In the following we make Monte Carlo simulations by 
considering 21 values of Ro from 0.1 to 10, sampled evenly 
on a logarithmic scale. We use the observed optical lumi- 
nosity as lo and keep redshift unchanged. Then lo, Ix and 
fluxes can be determined using Equations (3)-(5). The X- 
ray flux limit is determined as follows. Quasars in this sam- 
ple are from different observations and thus not uniformly 
sampled. All quasars with z ^ 3.5 are X-ray detected. In 
the range of « > 3.5, six quasars are not X-ray detected. 
Five of them, i.e. SDSS 1737+5828, PSS 1435-^3057, SDSS 
1532-0039, PSS 1506-^5220 and PSS 2344-^0342 were ob- 
served by Chandra with exposure times from 2.61-5.1 ks; 
SDSS 0338-1-0021 was observed by XMM-Newton with ex- 
posure time 5.49 ks. PSS 1506+5220 has one count in 0.5-2 
koV, and all the other flvo have zero counts. All the other 
detected sources have counts larger than 1. The average rest- 
frame /2 koV flux for one count in 0.5-2 keV of the six quasars 
is 0.8 X 10^'^^ erg cm~^ s~^ Hz~^. So we simply put a flux 
limit of /2 kev, upper = 1.5 X 0.8 X 10~^^ erg cm~^ s~^ Hz~^. 
We assume that simulated quasars with z > 3.5 and /2 kev 
less than /2 kev, upper will not be detected: quasars with 
0^/2 keV < 0.5/2 kev, upper wiU be assigned zero count, and 

quasars with O.5/2 keV, upper < /2 keV < 1.5/2 koV, upper wfll 

be assigned 1 count. Then the upper limits of non-detected 
quasars are at the 95% confidence level and will be calcu- 
lated according to Kraft et al. (1991), assuming one count 
corresponds to a flux of 0.8 x 10~^^ erg cm~^ s~^. 

Then wo simulate 100 samples for each of 21 different 
Ra values. For each R„, we compute the average slopes and 
Spearman correlation coefficients psp from the 100 simulated 
samples and display the results in Figures 2 and 3. Parame- 
ters for simulated samples are calculated in the same way as 
for panels (a) and (c) in Figure 1. As discussed in Strateva 
et al. (2005), the measurement errors and variability effects 
bring ax ~ 0.23 and ao ~ 0.17, and y/a'i + cr'i ~ 0.29. 
The observed dispersion is 0.4 for the M06 sample and 
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Figure 2. Spearman rank correlation coefficients psp as a mea- 
sure of the aox ^lo and aox ^Ix correlation for Miyaji's sample 
and simulated samples. Open squares: aox~lx in simulated sam- 
ples; solid circles: aox — lo in simulated samples; upper dashed 
line: aox — Ix of true data; lower dashed line: aox — lo of true 
data; solid line indicates no correlation. 
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Figure 3. Slopes of the aox — lo and aox — Ix correlation for 
Miyaji's sample and simulated samples. Symbols and lines are as 
in Figure 2. 

around 0.35 — 0.4 for other samples. The extra dispersions 
could be assigned to either Ox or (Jo as unknown disper- 
sions, so the possible range of R^ should be 0.5 ~ 1.4, i.e. 
log(J?a) ~ -0.3 - 0.15. In the range of log(i?<,) ~ -0.2 - 0, 
Psp and slopes in simulated samples for both aox — Ix and 
aox — lo correlations are consistent with observed values 
within 1.5(T, as shown in Figures 2 and 3. 

From the above two analyses, we conclude that a con- 
stant aox which does not depend on optical luminosity is 
still consistent with data in this high-z sample. 

2.2 Determining factor for an artificial correlation 
caused by luminosity dispersions 

Yuan et al. (1998) pointed out that a dispersion larger for 
the optical luminosity than for the X-ray luminosity tends 
to result in apparent, yet artificial correlation of aox — lo- 
To quantitatively examine the effect of dispersions on the 
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artificial correlation, wc make a simple analytic calculation 
as follows. 

As shown in Figure 4, assuming = Zo + const., disper- 
sions in and lo are and cto, respectively, and the span 
range of lo is Alo- Then locations of AGNs in the aox — lo 
plane with average Ix and la range of lo are indicated by 
the dotted lines in the lower panel, where an apparent cor- 
relation appears. We only consider the simplest situation: 

(i) AGNs are evenly distributed along lo (A concentration 
around a central lo is equivalent to a smaller Alo); 

(ii) AGNs arc only distributed 0.3838/(To away from aox, 
i.e. aox = aox ±0.3838/cto, where / is an unknown positive 
coefficient to be determined. We will discuss the value of / 
later. 

Then we fit the observed AGNs in the aox^lo plane, assum- 
ing a least chi-squared fitting procedure with same weight 
in lo and ao (i.e. aox/0.3838) as 

x' = [^((/oW-/c,(i))'+(aox(i)-aox(i))V0-3838')]/<T^, 

i 

. (6) 

where the lo(i) and aox{i) arc observed values, lo{i) and 
ceox (i) arc the values in the fitting line (indicated by a long 
solid line with a negative slope in the lower panel of Figure 4) 
with the least (/„ (i) - + (aox («) - aox (i))^/0.3838^ , 

and ffc is the typical constant error of lo- Then 



2 

X = 



k A/^ 1 

5(^(Y+3AU/.„f(l + -)^)/.5, (7) 



where k is the slope of the fitting line. The best fit slope k 
can be derived by solving 



dk ' 



The solution is 



All 



24/2 <t2 



+ 1 



(8) 



(9) 



(10) 



The deviation of the approximation in Equation (10) 
from Equation (9) is less than 20% when "^"^1^^^° < 1. Con- 
sidering two sub-samples with the same weight and differ- 
ent slopes fci and k2 respectively, the slope of combined 
sample including all points in each sub-samples would be 
fee = tan((arctan(A:i) -I- arctan(A;2))/2) ~ (fci + fc2)/2, where 
the last ~ is valid only if fei ^ 1 and A:2 -C 1. Therefore, 
we can derive the slope of the combined sample considering 
difTerent / values with different weights by an integration 



roo 

k = / kfpU)df, 
Jo 



(11) 



where p{f) is the probability of /. 

Then we take the distribution of / into account. If ax = 
0, and (Jo follows Gaussian distribution, the probability of / 
would be 



P(/) = 



0. 



(12) 
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Figure 4. Schematic sketches for the Ix — lo and aox — lo re- 
lationsliip. Here Ix = lo + aox/0.3838 and aox = —1.4. See 
section 2.2 for details. 



Then 

k 



2 -4- c/ '^o n2 



12ai 



^^^-,/(^)^ + l. 



12ai' 



(13) 
(14) 



Equation (13) shows that the slope of the artificial 
aox — lo correlation is directly proportional to al/Al^. The 
significance of the correlation, which could be measured by 
Spearman correlation coefficient, is always positively corre- 
lated to the slope value in the artificial aox — lo correlation, 
as shown in Figures 2-3. Therefore, the result in Equation 
(13) also means the significance of the artificial aox — lo 
correlation is proportional to a1/ All. 

For the M06 sample, Alo ^ 1.7 and CTo ~ 0.40 when 
ax = 0, then k ~ —0.3. Such estimation of an artificial 
slope using Equations (13) or (14) is qualitatively consistent 
with the slope in the e = 1 simulations, where the value is 
k ~ —0.18 as shown in Figure 3. A reason for the discrepancy 
is that the absolute slope value is not <IC 1, therefore ap- 
proximations used in our estimation are deviated from true 
values. Another possible reason, i.e. different fitting proce- 
dure used in simulations and our estimation, would more 
or less contribute to the lower absolute slope values in sim- 
ulations. In simulations, a linear regression method, which 
only takes residuals of the dependent variable into account, 
is used, thus always leading to a lower absolute slope value 
than methods considering residuals in both variables as used 
in Equation (7), as shown in Panel (b) in Figure 1. 

We now discuss consequences of non-zero ax- When 
ax > 0, the distribution of / will be extended. As shown 
in Figure 4, point B becomes a distribution in the range of 
CD within la. When calculating the x'^ of a linear regres- 
sion with slope —k, the contribution of the broadening in B 
is equivalent to an extension along the lo axis. For example, 
as shown in Figure 4, the contribution of point C is equiva- 
lent to point C which has the same aox as B but smaller lo. 
Therefore, a non-zero ax tends to smooth the distribution of 
lo and extend its range. When ax/k is comparable with lo, it 
will extend the range of lo significantly. To show this effect, 
we do another simulation. Based on the M06 sample with a 
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Figure 5. Spearman rank correlation coefficients psp and best- 
fit slopes of tlic ciox — lo correlation for simulated samples for 
various Ux, with constant co = 0.3 and aox = —1.729. lo and 
redshift distributions are from Miyaji et al. (2006). 

given (Jo = 0.3, we calculate the Spearman correlation coef- 
ficients and slope for aox — lo correlation with different , 
while other conditions are set to be the same as simulations 
in section 2.1. As shown in Figure 5, when log{ax/(Jo) > 0, 
i.e. > 0.3 and a^/k > 1.7 ~ Alo, both the Spearman cor- 
relation coefficient and the slope tend to move toward zero 
when ax increases. When \og{ax/cro) > 0.5, i.e. ffx > 0.9 
and ax/k > 5.3 ^ 3Alo, there is no correlation in aox — h 
within la. However, as discussed in Section 2.1, even if all 
extra dispersion in aox comes from Ix, ax is unlikely to 
exceed 0.4 and \og{ax/ao) is unlikely to exceed 0.3. 

The effects of ao on the relationship of aox — Ix is 
similar with the effects of ax on the relationship of aox — lo- 
Thus we do not need to repeat the above analysis. Moreover, 
a similar effect as presented here for the aox — lo correlation 
would also affect any correlation with a dependent variable 
B, which is not directly observed but derived from B oc A'^, 
where A is the independent variable, such as the Baldwin 
effect, which has also been pointed out by Yuan et al. (1998) 

In summary, the significance of artificial correlation in 
aox — lo is approximately proportional to ctq/AZo, and de- 
creases when ax increases and becomes comparable with 
kAlo, where k is the absolute value of the artificial slope 
and Alo is the range lo span. 



3 DATA ANALYSIS OF A SUB-SAMPLE OF 
AGNS FROM STEFFEN ET AL. (2006) 

Another problem in the study of aox — lo relationship is the 
degeneracy between redshift and luminosity in flux-limited 
samples. S06 used a much larger sample than M06 with 
Alo ~ 5, which suppresses the false slope artifacts discussed 
in section 2. The observed change in aox across this larger 
baseline Alo in their sample is sufficiently large that aox 
must depend on luminosity, or redshift, or both. To distin- 
guish between luminosity dependence and redshift depen- 
dence, first, S06 performed partial correlation analysis us- 
ing Kendall's generalized partial t to quantitatively show 
the correlation significance of aox — lo and aox — z. They 
found a 13.6cr correlation of aox — lo when controlling z, 



and a 1.3o" correlation of aox — z when controlling lo- How- 
ever, KoUoy et al. (2007) has showed that interpretation of 
Kendall's r is problematic, and Kendall's r for the aox — z 
correlation is not necessarily expected to be non-zero when 
aox is correlated with z. Based on simulations, Kollcy et 
al. (2007) pointed out that the lack of evidence for a sig- 
nificant correlation between aox and z based on Kendall's 
T in StefFen et al. (2006) may be the result of am incor- 
rect assumption about the distribution of t under the null 
hypothesis. In spite of most previous studies, KoUey ct al. 
(2007) found that aox is correlated with both lo and z. 
Moreover, in the partial correlation analysis of aox — lo, 
consequences of luminosity dispersions, as discussed in sec- 
tion 2, were not taken into account. To show this effect, we 
select AGNs in three redshift bins, with each bin contain- 
ing 38 sources, to control the redshift. Then we examine the 
aox — lo and aox — Ix relations in each bin, as shown in 
Figures 6 and 7. Similar to Figure 1, in each redshift bin, 
aox is anti-correlated with lo, but positively correlated with 
Ix, which is caused by luminosity dispersions. Because dis- 
persions always strengthen the anti-correlation of aox — lo, 
dependence of aox on lo might be biased toward higher sig- 
nificance by luminosity dispersions, whereas the dependence 
on z does not suffer such bias. 

Second, S06 compared aox residuals as a linear func- 
tion of lo and z. As shown in their Figure 8, there are sys- 
tematic residuals of aox — aox (z), which indicate that aox 
cannot be linearly dependent on redshift alone. However, 
spectral index might depend on redshift in a non-linear form, 
as shown in Figure 12 of Strateva et al. (2005). Moreover, it 
is also possible that aox depends on both lo and z, as shown 
in Kelley et al. (2007). Using different parameteric models 
for the redshift and optical luminosity dependencies, Kelley 
et al. (2007) found the model that is best supported by their 
data has a linear dependence of aox on cosmic time, and a 
quadratic dependence of aox on lo (the definition of aox 
in Kelley et al. (2007) is different from our definition with 
an opposite sign). Since lo and z are coupled together in 
flux limited samples, different parameteric models will lead 
to different results and their best model results depend on 
the form of the models. In summary, dependence of aox 
on z, though with lower significance in partial correlation 
analysis, could not be excluded in S06 sample. 

To avoid possible bias from aox — z correlation, here 
we examine a sub-sample from S06 containing 187 AGNs, as 
shown in the dotted-line box in Figure 3 in S06 (Steffen et al. 
2006; Strateva et al. 2005; Vignah et al. 2005; Shemmer et 
al. 2005; Kelly et al. 2005), which more completely fills the 
redshift and optical luminosity plane and thus is less affected 
by such bias. We refer this sample as the 'low-z sub-sample'. 
We perform linear regressions in the same procedure as de- 
scribed in section 2.1. Figure 8 presents our results for this 
sample. Similar to the M06 sample, the low-z sub-sample 
show conflicting correlations due to dispersions in luminos- 
ity; aox = — O.I6Z0 -I- const as shown in solid line in Panel 
(a), but aox = 0-lllx + const as shown in solid line Panel 
(c). Therefore the same data produce two totally different 
results: Ix = 0.58Zo + const or Ix = 1.40Zo + const, i.e. e < 1 
or e > 1 if Lx oc Lfj. As shown in Panel (b), the slope of the 
Ix — lo relation also depends on which luminosity is used as 
the dependent variable. When treating lx as the dependent 
variable, the slope is 0.59 ± 0.06 (the flatter dashed line), 
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Figure 6. aox dependece on the rest frame 2500 A monochro- 
matic luminosity for AGNs in three redshift bins in S06 sample: 

z < 0.1 (stars), 1.45 < z < 1.77 (circles) and 4 < z < 4.6 
(downward-pointing triangles). X-ray detected AGNs are repre- 
sented using filled symbols while upper limits arc represented us- 
ing open symbols. The solid line, dashed line and dash-dotted 
line indicate linear regression results from the EM Algorithm in 
ASURV for z ^ 0.1, 1.45 < 2 < 1.77 and 4 < 2 < 4.6 samples, 
respectively. 
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Figure 7. apx dependece on the rest frame 2 keV monochro- 
matic luminosity for AGNs in three redshift bins in S06 sample; 
z ^ 0.1 (stars), 1.45 < 2 < 1.77 (circles) and 4 < 2 < 4.6 
(downward-pointing triangles). X-ray detected AGNs are repre- 
sented using filled symbols while upper limits are represented us- 
ing open symbols. The solid line, dashed line and dash-dotted line 
indicate linear regression results using X-ray detected sources for 
2 ^ 0.1, 1.45 < 2 < 1.77 and 4 < 2 < 4.6 samples, respectively. 



while it changes dramatically to 1.45 ± 0.12 (the steeper 
dashed line) when treating lo as the dependent variable. Us- 
ing ILS bisector, we find the slope to be 0.93 ± 0.08 (solid 
line), which is consistent with e = 1 (dotted line). In panel 
(c), the fit looks different from the trend by eye (bisector 
fit). As discussed in section 2, when using traditional ordi- 
nary least-squaxed method which minimizes the residuals of 
the dependent variable, the fit tends to be flatter than the 
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Figure 8. 187 low-redshift AGNs from Steffen et al. (2006). Cir- 
cles indicate X-ray detected AGNs, while arrows indicate up- 
per limits. Panel (a) apx dependence on the rest frame 2500 
A monochromatic luminosity. The solid line indicates linear re- 
gression results from the EM Algorithm in ASURV. Panel (b) 
rest frame 2 keV monochromatic luminosity against 2500 A one. 
The solid line indicates the ILS bisector result, and dashed lines 
indicate ILS(Y|X) and ILS(X|Y) results from the EM Algorithm 
in ASURV, respectively. The e = 1 relation is shown by dotted 
line for comparison. Panel (c) qqx dependence on the rest frame 
2 keV monochromatic luminosity. The solid line indicates linear 
regression result using X-ray detected quasars. Spearman corre- 
lation coefficients and slopes of fitting lines are indicated in the 
upper right panel. 



bisector one where residuals of both dependent and indepen- 
dent variables are taken into account. When data points are 
concentrated at the center, as in Figure 8(c), inconsistence 
of the two fits becomes large. 

Interestingly, but not surprisingly, results of ILS (aox | 
lo) are consistent with results of ILS(?2! j lo), and results 
of ILS(aox I Ix) are consistent with results of ILS(Zo | Ix), 
as shown in Table 1 and Figures 1 and 8. The reason is 
that artificial correlations seen in aox — lo relation, i.e. 
0.3838(^1: — lo) — lo relation, corresponds to least-squares 
method which minimizes the residuals of Ix, which is the 
same as in Ix — lo relation, and thus leads to a slope less than 
the true one where residuals of both variables are considered. 
Such results indicate that for Ix — lo relation, regression re- 
sults based on the traditional ordinary least-squares method 
which minimizes residuals of the dependent variable sufi^er 
from effects caused by luminosity dispersion and are thus 
not reliable. Instead, weighting both lo and Ix in the regres- 
sion, i.e. ILS bisector, would be a more robust methodology. 
Moreover, the degree to which the slopes of ILS(aox | Ix) 
and ILS(aox | lo) are inconsistent indicates the degree of 
artificial correlation discussed in section 2. When the slopes 
converge to the same value, the artificial correlation would 
be suppressed. 

We conclude for this sample that a constant aox which 
does not depend on luminosity is also consistent with data. 
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Table 1. Slopes of Ix — lo relation in M06 and low-z sub-sample 
derived from different regression methods. 



Regression method 


M06 


low-z sub-sample 


ILS(aox 1 lo) 


0.53 ± 0.13 


0.58 ± 0.05 


ll^^x 1 lo) 


0.54 ± 0.14 


0.59 ± 0.06 


ILS(aox 1 «x) 


2.08 ± 0.48 


1.40 ± 0.10 


ILS(/o 1 Ix) 


2.08 d= 0.45 


1.45 ± 0.12 



4 COMPARISON OF OPTICAL AND X-RAY 
LUMINOSITY FUNCTIONS 

For a complete sample of broad line AGN including op- 
tical and X-ray observations (detections or upper limits), 

slopes in optical and X-ray quasar luminosity functions 
(LPs) should be the same after correct transformations. 
When converting optical luminosity to X-ray luminosity, dif- 
ferent aox models would lead to different X-ray LF shapes 
in the optical frame. Thus we can test whether a particular 
aox model is correct by comparing the two LFs. Wo use 
the optical quasar LF from Richards et al. (2006) and AGN 
hard X-ray LF from Barger et al. (2005). 

Wc investigate the following two aox models: 

(i) a constant aox, Ix =lo — 3.92; 

(ii) aox from Steffen et al. (2006), Ix = 0.72lo + 4.53; 

We follow Hopkins, Richards & Hernquist (2007) to cal- 
culate the binned LFs. For each aox model, an overall nor- 
malization factor is applied in the X-ray LF to get the mini- 
mum x^, which means that we are comparing just the slopes 
of LFs. Results are presented in Figure 9. A constant aox 
which does not depend on luminosity (left panels) is consis- 
tent with data, and is more preferred than the aox models 
given by S06 (right panels). 

However, from this comparison we cannot reach a strong 
conclusion that luminosity dependent aox is excluded com- 
pletely, for three reasons as follows. First, there are very few 
data points hero in tho X-ray LF. Second, as pointed out by 
Richards et al. (2005) , such comparison is not strictly quan- 
titative since X-ray selected samples and optically selected 
samples are not identical. Moreover, the bright-end slopes 
from different X-ray samples are different (Barger et al. 
2005; Ucda ct al. 2003; Hasingor, Miyaji, & Schmidt 2005). 
Hopkins, Richards & Hernquist (2007) combined a large set 
of LF measurements and took obscuration and scattering 
into account, and in their analysis the luminosity functions 
can be reconciled reasonably well with the aox model in 
S06. However, since the constraints on the present bright- 
end X-ray LFs are poor, the fact that the LFs could be 
reconciled reasonably well with S06 in their work probably 
just reflects the large X-ray error bars. 



5 SELECTION EFFECTS IN FLUX LIMITED 
SAMPLES: OPTICALLY SELECTED 
SAMPLES VS X-RAY SELECTED SAMPLES 

In this section, wc discuss the selection effects in flux limited 
samples. For a given AGN luminosity function, assuming the 
observed optical and X-ray luminosities are the intrinsic val- 
ues modifled by dispersions which might be caused by vari- 
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Figure 9. Comparison of the SDSS DR3 optical quasar luminos- 
ity function from Richards ct al. (2006) (lines with errorbars) with 
AGN hard X-ray luminosity function from Barger et al. (2005) 
(crosses with errorbars). Left panels show the LFs with e = 1 
(luminosity independent spectral index). Right panels show the 
results with e = 0.72 (spectral index is anti-correlatcd with lu- 
minosity). Different rows show results at different redshift range. 
For each panel, an overall normalization factor is applied in X-ray 
LF to get the minimum x^- 

abilities or observational errors, there are three possibilities 
in a flux limited sample: 

a) lower fraction of more luminous AGNs axe missed; 

b) higher fraction of more luminous AGNs are missed; 

c) same fractions of more luminous and fainter AGNs 
are missed, so the relationship between Ix ~ lo remains un- 
changed. 

Assuming the slope of Ix — lo relation (without fur- 
ther description, slopc=e in Ix = e x h + const throughout 
this section) is unity, and at a certain redshift the rmmbor 
density of AGN decreases with luminosity. Figure 10 shows 
schematic sketches for the first two cases in optically selected 
AGN samples. The upper panels are for case (a), where the 
density contour lines in luminosity functions in the lo — z 
plane (each lino corresponds to a given constant AGN num- 
ber density as a function of redshift) are steeper than the 
flux limits, hence lower fraction of more optical luminous 
AGNs are missed, and then the slope is biased toward more 
than unity. The bottom panels are for case (b), where the 
density contour lines in luminosity functions in the lo — z 
plane axe flatter than the flux limits, hence higher fraction 
of more optical luminous AGNs are missed, and then the 
slope is biased toward less than unity. The loft panels show 
flux limits (dashed lines) in lo — z plane, compared with the 
density contour lines in luminosity functions (solid lines). 
The left panels show observed AGNs (solid circles) which 
are above the flux limits (deished lines), and missed AGNs 
(open circles) which arc below the flux limits, in lo — lx plane. 
Slope= 1 are indicated by solid lines. 

The density contour lines in real luminosity functions 
are much more complicated than shown in Figure 10, and the 
real slopes of density contour lines depend on redshifts and 
luminosities. Moreover, since the fraction of missed AGNs 
depends on the dispersions of luminosities around the linear 
relationship of Ix — lo, the biases also depends on the disper- 
sions. To test whether the slope of Ix — lo relation could be 
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Figure 10. Schematic sl^etches for selection effects in optically 
selected AGN samples. The upper panels are for case (a) and 
the bottom panels arc for case (b). The left panels show flux 
limits (dashed lines) valo — z plane, compared with the density 
contour lines in luminosity functions (solid lines). The left pan- 
els show observed AGNs (solid circles) which arc above the flux 
limits (dashed lines), and missed AGNs (open circles) which are 
below the flux limits, in lo — Ix plane. Slope is defined as in Ix — 
slope xio + const. Slope= 1 are indicated by solid lines. Note that 
in the upper-right panel slope> 1 for detected AGNs, and in the 
lower-right panel slopc< 1 for detected AGNs, since Ix is the x- 
axis and lo is the y-axis. These sketches, much simplified and only 
qualitatively correct, are shown for illustration only. 

biased by flux limits in realistic optically selected AGN sam- 
ples, we carry out Monte Carlo simulations. We use optical 
analytical luminosity function from Richards et al. (2005) 
for a < 3 AGNs, and Richards et al. (2006) for « > 3 AGNs. 
We simulate three optically selected sub-samples, in order 
to mimic the SDSS, COMBO-17 and high-z samples in S06: 

1) a shallow z < 3 sub-sample with rUg < 19 containing 
about 155 AGNs, which is similar to the SDSS sample in 
S06; 

2) a deeper z < 3 sub-sample with rUg < 21 containing 
about 52 AGNs, which is similar to the COMBO-17 sample 
in S06; 

3) a 4 < 2 < 6 sub-sample with nig < 20 containing 
about 55 AGNs, which is similar to the high-z sample in 
S06. 

We do not try to mimic the nearby Seyfert 1 and BQS 
samples in S06 which are located in z < 0.4, since the LF in 
this redshift range has larger errors due to smaller volume. 
For the two z < 3 sub-samples, a detection efficiency factor 
is taken from Figure 6 in Richards et al. (2006). For the z > 4 
sub-sample, a constant detection efficiency is used according 
to Richards et al. (2006) . To show the effects of luminosity 
dispersions in optically selected flux limited samples, assum- 
ing the slope of Ix — lo equals unity with dispersions, 1000 
simulations, each containing the above three sub-samples, 
are carried out for each of the following five dispersion mod- 
els: 

(i) ao = 0.42, ax = 0; 

(ii) ao = ax = 0.3 for 2; < 3 AGNs, cTq = fx = 0.5 for 
2 > 4 AGNs; 

(iii) ao = ax = 0.3 

(iv) ao = ax = 0.5 



(v) ao = 0,ax= 0.42. 

where ao and ax are defined as in section 2.1. Models (i) and 
(v) are extreme cases, in order to show the bias direction 

when ao or ax is dominating. Models (ii) ^ (iv) can show 
the effects of dispersion, and dispersion evolution in cosmic 
time. 

The Monte Carlo analysis was performed by generating 
a combined sample containing the above three sub-samples 
for each of the five dispersion models as follows: first, the 
redshift and luminosity ranges are divided into grids with 
5z = 0.1 and 5Mg = 0.3, then the detection probability in 
a given grid is proportional to dV{z) x ^{z, Mg) x rj, where 
dV{z) is the volume element in comoving space, ^{z,Mg) 
is the luminosity function, and rj is the detection efficiency. 
In a given grid, AGNs axe randomly produced following an 
uniform distribution, and the total number of AGNs in the 
grid is determined by a poisson process with expectation 
N{z, Mg) = dV{z) X ^{z, Mg) xr/xC, where C is a constant 
for a given sub-sample with given luminosity dispersions, 
adjusted to make the average number of detected AGNs to 
be 155, 52 and 55 for the three sub-samples, respectively. 
The conversion from A4g to the intrinsic optical luminos- 
ity lo follows Richards et al. (2005), and the intrinsic X-ray 
luminosity Ix = lo — 4. Second, luminosity dispersions are 
applied, where the observed lo and Ix are drawn from the 
Gaussian distribution around lo and Ix with given disper- 
sions ao and ax, respectively. Third, flux limits are applied, 
AGNs with nig > flux limit are detected. The above pro- 
cedure select about 262 AGNs for each dispersion model, 
where about 155 in sub-sample 1), 52 in sub-sample 2) and 
55 in sub-sample 3). Then the above procedure is repeated 
1000 times to get 1000 independent samples. 

The slope in each simulated sample is calculated us- 
ing FITEXY (Press et al. 1992) assuming the same error 
in lo and Ix- The distributions of slopes in simulated opti- 
cally selected flux limited samples are shown in Figure 11. 
The relative probabilities are normalized with peak val- 
ues equal unity. From left to right are the five dispersion 
models (i) to (v) respectively. The mean values ± stan- 
dard deviations of slopes in simulations of the five models 
are: 0.76 ± 0.05, 0.89 ± 0.06, 0.97 ± 0.05, 1.04 ± 0.11, and 
1.23 ±0.06, respectively. In models (iii) and (iv), ao = ax = 
0.3 and 0.5, flux limited simulations are consistent with the 
assumed slope= 1, i.e. the slopes are not biased. However, 
when (Jo 7^ ax or ao, ax change in cosmic time, slopes in flux 
limited simulations might be biased, as shown in models (i), 
(ii) and (v). Slopes in model (i), i.e. 0.76 ± 0.05, is consis- 
tent with the slope in S06, i.e. 0.72 ± 0.01. However, since 
(To = 0.42 and ax = is an unrealistic extreme case, this 
does not mean that the non-unity slope of Ix — lo relation 
is totally caused by such selection effects. Note that slope 
values in Figures 11 depend on the selection method, i.e. 
flux limits and number of AGNs in each sub-sample, hence 
another different optical samples will have different results. 

We do not simulate X-ray selected flux limited samples 
for two reasons. First, if the slope of lo — Ix equals unity. 
X-ray AGNs are identical to optical AGNs, and the X-ray 
LF will be the same as optical LF with Ix = I0+ constant, 
as suggested in section 4. Therefore, the results here can 
be applied to X-ray selected sample with similar flux limits 
after switching ao and ax- Second, X-ray LFs have larger 
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Figure 11. The distributions of slopes in simulated optically 
selected flux limited samples, where Ix =slopexio+constant. The 
relative probabilities are normalized with peak values equal unity. 
^From left to right arc the five dispersion models (i) to (v) listed 
in the text, respectively. 



errors than optical LFs due to smaller samples. Therefore, 
it is our purpose to just point out the fact that the slope will 
be biased in flux limited X-ray samples, rather than focusing 
on a particular X-ray selected sample. 

While a number of previous studies of optical selected 
AGNs have reported that aox is anti-correlated with lumi- 
nosity (e.g. Strateva et al. 2005; Steffen et al. 2006), whether 
aox depends on luminosity in X-ray selected AGN samples 
remains unknown. Hasingor (2004) found no aox depen- 
dence on either luminosity or redshift in soft X-ray selected 
samples. Prank et al. (2007) found in their Chandra Deep 
Field-North sample = (0.808±0.047)Zo -fconst. They also 
found the slope decreases when only brighter sources arc in- 
cluded, and the slope increases when only fainter sources 
are included. It is possible that the slope might be biased 
in this flux limited sample and the magnitude of biases are 
different when using different flux limit. As shown in Fig- 
ure 11, if (To > cTx, an optically selected sample like in S06 
will be biased toward a flatter slope. Moreover, if the X-ray 
LF is similar to optical LF and the X-ray sample is consisted 
of AGNs with similar flux limits, the X-ray sample will be 
biased toward a steeper slope. Therefore, even if optically 
selected AGNs and X-ray selected AGNs arc identical, the 
slope in the optically selected sample will be flatter than 
the slope in the X-ray selected sample, which can properly 
explain the observed discrepancy. 

In summary, selection effects in flux limited samples 
might bias the lo — Ix relation and cause discrepancy in the 
lo — Ix relation in optically selected samples and X-ray se- 
lected samples, especially when cTq =^ ffx or cto, cTx change 
in cosmic time. The magnitude of the bias and discrepancy 
depend on the luminosity function, flux limits of the sam- 
ple, and dispersions in optical and X-ray luminosities. Note 
that even if such selection effects do bias the slope of the 
lo — Ix relation toward the observed discrepancy between 
optical and X-ray samples, it is not necessarily the only rea- 
son. It is possible that optically selected samples and X-ray 
selected samples are consisted of different AGNs, so slopes 
of the lo — Ix relation in optically selected samples will be 



different from slopes in X-ray selected ones. As discussed in 
Brusa et al. (2007), about 40% of the X-ray selected AGNs in 
their COSMOS sample would have not been easily selected 
as AGN candidates on the basis of purely optical criteria, 
either because similar colors to dwarf stars or field galaxies, 
or because they are not point like sources in morphologi- 
cal classification. Moreover, optically selected AGNs and X- 
ray selected AGNs might be typically in different evolution 
stages and thus are not identical (Shen et al. 2007). 



6 DISCUSSION AND CONCLUSIONS 

In summary, we have investigated the correlation between 
the spectral index aox and optical/X-ray luminosities in 
AGNs by means of linear regressions, Monte Carlo simular- 
tions, simplified analytic estimations and comparison of X- 
ray and optical luminosity functions. We have reached five 
conclusions: 

1. The dependence of aox on optical luminosity found 
in Miyaji et al. (2006) may not be an underlying physical 
property. It remains unknown whether e < 1 or e > 1 if 
Lx oc Lq in this high-z sample. 

2. The luminosity dependence can be artificially gener- 
ated very easily by luminosity dispersions. The significance 
of artificial correlation in aox — lo is approximately propor- 
tional to ao/Al'o, where ctq is the optical luminosity disper- 
sion and Alo is the range that lo spans, and decreases when 
ax increases and becomes comparable with kAlo, where k 
is the absolute value of the artificial slope. This effect also 
affects the Baldwin effect. Instead of regressions only weight- 
ing one variable, weighting both lo and lx, i.e. ILS bisector, 
in the regression would be a more robust methodology to 
avoid such bias. 

3. In a more complete low-z sub-sample from Steffen et 
al. (2006), aox must depend on luminosity, or redshift, or 
both. However, a luminosity independent aox is still consis- 
tent with data. Redshift dependencies cannot be ruled out 
and may be large, but somewhat hidden because of lumi- 
nosity dispersions, which generate artificial luminosity cor- 
relations in each redshift bin. 

4. In the comparison of X-ray (Barger et al. 2005) and 
optical quasar (Richards et al. 2006) LFs, a luminosity in- 
dependent aox is consistent with data, and more preferred 
than the luminosity dependent aox model given by S06. 

5. Selection effects in flux limited samples might bias 
the lo — lx relation and cause discrepancy in the lo — lx rela- 
tion in optically selected sample and X-ray selected sample, 
especially when ao 7^ (Jx or tJo, (Jx change in cosmic time. 
The magnitude of the bias depends on the luminosity func- 
tion, flux limits of the sample, and dispersions in optical and 
X-ray luminosities. 

It therefore remains inconclusive whether the anti- 
correlation between AGN spectral index and optical lumi- 
nosity is true. Even if aox does depend on optical luminos- 
ity, the currently adopted slope value might be biased and 
deviate from the intrinsic value. To correctly establish a de- 
pendence of aox of AGNs on their luminosity, a larger and 
more complete sample, such as from multi-wavclcngtli sur- 
veys, is needed and consequences of luminosity dispersions 
and selection effects in flux limited samples must be taken 
into account properly. 
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